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The monodromy transform and corresponding integral equation method described here give rise to 
a general systematic approach for solving integrable reductions of field equations for gravity coupled 
bosonic dynamics in string gravity and supergravity as well as for pure vacuum gravity in four 
and higher dimensions. For physically different types of fields in space-times of D > 4 dimensions 
with d = D — 2 commuting isometries - stationary fields with spatial symmetries, interacting waves 
or evolution of partially inhomogeneous cosmological models, the string gravity equations govern 
the dynamics of interacting gravitational, dilaton, antisymmetric tensor and any number n > of 
Abelian vector gauge fields (all depending only on two coordinates). The equivalent spectral problem 
constructed earlier allows to parameterize the entire infinite-dimensional space of (normalized) local 
solutions of these equations by two pairs of arbitrary coordinate-independent holomorphic d x d- 
and d x n- matrix functions {u±(w), v±(w)} of a spectral parameter w which constitute a complete 
set of monodromy data for normalized fundamental solution of this spectral problem. The "direct" 
and "inverse" problems of such monodromy transform — calculating the monodromy data for any 
local solution and constructing the field configurations for any chosen monodromy data always 
admit unique solutions. We construct the linear singular integral equations which solve this inverse 
problem. For any rational and analytically matched (i.e. u+(w) = u_(w) and v+(iu) = v_(w)) 
monodromy data the solution of these integral equations and corresponding solution for string 
gravity equations can be found explicitly. Simple reductions of the space of monodromy data leads 
to the similar constructions for solving of other integrable symmetry reduced gravity models, e.g. 
5D minimal supergravity or vacuum gravity in D > 4 dimensions. 



PACS numbers: 04.20.Jb, 04.50.-h, 04.65.+e, 05.45.Yv 
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Introduction 

Motivated by various physical reasons, numerous mod- 
ern studies of strong gravitational fields in different grav- 
ity models - pure gravity as well as the bosonic sectors 
of string gravity and supergravity in four and higher di- 
mensions, brought us many interesting discoveries con- 
cerning the existence and properties of a large variety of 
higher-dimensional space-time structures - black holes, 
black rings, black lens, etc. (see, for example, [l|), which 
distinguish these models in higher dimensions from those 
in 4D space-times. It can be expected that the same may 
concern also physically different types of fields, e.g. in- 
teracting waves and cosmological solutions. In all cases, 
an explicit construction of appropriate solutions can play 
even more crucial role than it was in four dimensions, 
where our physical intuition is better adapted to. 

The construction of solutions for different gravity mod- 
els in four and higher dimensions needs a generalization 
of various methods developed earlier for 4D case to make 
these available for solving the field equations in higher di- 
mensions and to include various non-gravitational fields 
which play important roles in these theories. For these 



purposes, in particular, it is very important to find the 
cases in which the symmetry reduced dynamical equa- 
tions occur to be completely integrable. Such integrable 
reductions may arise if space-times of D dimensions ad- 
mit D — 2 commuting isometries and therefore, all field 
components and potentials depend only on two space- 
time coordinates. Now we have a large experience in the 
development and application of various solution generat- 
ing methods for Einstein's field equations in AD space- 
times with two commuting isometries (see, e.g., the sur- 
vey [H and the references therein). In AD case, the al- 
ready developed methods allow, in particular, to con- 
struct explicitly the hierarchies of solutions with arbi- 
trary large number of free parameters. These are the 
./V-soliton solutions generated on (or, in another words, 
interacting with) arbitrary chosen backgrounds such as 
Belinski and Zakharov vacuum solitons |3| and Einstein 
- Maxwell solitons |4j, as well as some classes of non- 
soliton solutions 0-0] • Besides that, for these equations 
one can solve the (effectively - two-dimensional) bound- 
ary value problems for stationary fields (see [8[ and the 
references there) and characteristic initial value problems 
However, a generalization of these methods for 
higher dimensional gravity models is not trivial and it 
may need to overcome some specific difficulties. 
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In the last two decades, many authors used the Belin- 
ski and Zakharov soliton generating transformations for 
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constructing asymptotically flat stationary axisymmet- 
ric solutions for pure vacuum 5D gravity. This "dress- 
ing" method for generating solitons on arbitrary (vac- 
uum) background can be generalized directly from AD 
to higher dimensions without any essential changes of 
the procedure. However for appropriate choice of soli- 
ton parameters and practical constructing of particular 
solutions with desired propeties for the fields of rotating 
objects in 5D space-times, the analysis of the so called 
"rod structure" of soliton solutions suggested in [l^, EH 
was found very helpful. (Note that many decades ago, 
the similar diagrams were used for classification of static 
vacuum solitons in AD case [l4j|.) 

A series of attempts was made also to apply the Belin- 
ski and Zakharov construction of the spectral problem to 
other gravity models in AD as well as to bosonic sector of 
5D minimal supergravity. The main problem which arose 
in these applications and which was not solved there is 
that the constructions of the Lax pairs (like the Belinski 
and Zakharov one) had not been supplied by appropriate 
conditions imposed on the fundamental solutions of the 
linear systems which would provide the equivalence of 
the entire spectral problems to the dynamical equations. 
In these cases, the application of Belinski and Zakharov 
soliton generating transformations was not successful be- 
cause the generating soliton solutions left the correspond- 
ing coset space (see the comments in (l5j). 

On the other hand, the spectral problem constructed 
in [l6j is equivalent to the dynamical equations which 
govern the bosonic sector of heterotic string effective the- 
ory in the space-time of D dimensions with d = D — 2 
commuting isometries. (These equations include as dy- 
namical variables the metric, dilaton field, antisymmetric 
potential for the 3-form gauge field and any number n of 
Abelian vector gauge fields, all depending only on two 
space-time coordinates.) In 5D case with n = 1, certain 
constraints imposed on the potentials reduce this system 
and corresponding equivalent spectral problem to those 
for the bosonic sector of 5D minimal supergravity. 

As it was mentioned in [l6j , the spectral problem con- 
structed there possess the monodromy preserving prop- 
erties, which provide the base for generalization for this 
case of the monodromy transform approach suggested by 
the author in [l7| for solving the symmetry reduced vac- 
uum Eisntcin equations, electrovacuum Eistcin - Maxwell 
and Einstein - Maxwell - Weyl field equations in General 
Relativity (see also [II EH)- 

In this paper we present some sketch of the monodromy 
transform construction for heterotic string gravity model 
in D dimensions and 51? minimal supergravity including 

• the solution of a "direct" problem - a parametriza- 
tion of the whole space of local solutions of 
symmetry reduced bosonic equations for heterotic 
string gravity by pairs of matrix valued coordinate- 
independent holomorphic functions of spectral pa- 
rameter {u^ d (w), v^ xn (w)} which arise as the 
monodromy data of the normalized fundamental 
solution of the associated spectral problem and 



which characterize uniquely any local solution; 

• the solution of the "inverse" problem - a construc- 
tion of the local solutions for any given monodromy 
data. For solution of this problem we construct a 
system of linear singular integral equations which 
solution for any given monodromy data always ex- 
ists and is unique, and it allows to calculate all field 
components in quadratures; 

• explicit construction of infinite hierarchies of multi- 
parametric families of exact solutions for a special 
class of "analytically matched" monodromy data 
such that u + (w) = \i-(w) = u(w) and v + (w) = 
v_(u>) = v(w). These conditions are equivalent 
to regularity of the local solution at the boundary 
which consists of stationary points of the isometry 
group (the "axis" of symmetry). The monodromy 
data functions u(w) and v(w) for this class are 
closely related to the values of the matrix Ernst po- 
tentials on the "axis" . If u(w) and v(w) are chosen 
as any rational functions of w, the corresponding 
solution can be found explicitly. 

Finally we note, that the monodromy data play the role 
of "coordinates" in the infinite-dimensional space of local 
solutions of symmetry reduced dynamical equations. If 
we impose certain restrictions on the choice of the mon- 
odromy data, the integral equation method developed 
here for solving the for bosonic field equations in het- 
erotic string theory can be used also for constructing so- 
lutions for pure vacuum gravity in D dimensions and 5D 
minimal supergravity (the details are expected to be the 
subject of some later publications). 

Massless bosonic dynamics in string gravity 

The massless bosonic part of string effective action in 
space-times with D > A dimensions in string frame is 

S =J e" $ !.R^ + V„W M i - ±H MNP H MNP 

2 p=l J 

J 1 ) 

where M,N,... = 1,2, ...,£) and p = 1, ...n; Gmn 
possesses the "most positive" Lorentz signature. Met- 
ric Gmn and dilaton field <!> are related to the metric 
Gmn and dilaton $ in the Einstein frame as 

Gmn = e 2 *G M N, $ = (£>- 2)$. (2) 

The components of a three-form H and two-forms 

are determined in terms of antisymmetric tensor field 

Bmn and Abelian gauge field potentials Am^ as 

n 

HmNP = 3(9[A/-B A r P ] - A^m^Fnp]^), 

p=l 

Fmn^ — 2d[M^A r ]^ p \ Bmn = —Bnm- 
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Space-time symmetry ansatz 

We consider the space-times with D > 4 dimensions 
which admit d = D — 2 commuting Killing vector fields. 
All field components and potentials are assumed to be 
functions of only two coordinates x 1 and x 2 , one of which 
can be time-like or both are space-like coordinates. We 
assume also the following structure of metric components 



G 



_ [9^ 



1,2 



G a b 

while the components of field potentials take the forms 



B 



MN 





B al 



•l.w"" ( A % ) N ' 



We choose x ,x so that g^ v takes a conformally flat form 



ei 
e 2 



ci = ±l 
(2 = ±1 



where fix 1 *) > and the sign symbols e\ and e 2 allow 
to consider various types of fields. The field equations 
imply that the function a(x , x ) > is "harmonic" one: 

,„ 2 



det 1 1 G Q 6 



fd^a = 0, 



-eiC2- 



where rf v is inverse to ?7 M „, and therefore, the function 
P(x^) can be defined as "harmonically" conjugated to a: 



d a /3 = eea v d v a, 



1 

-1 



Using the functions {a, ft), we construct a pair (£,r)) of 
real null coordinates in the hyperbolic case or complex 
conjugated to each other coordinates in the elliptic case: 



rj = P - ja, 



J 



1, e = 1 — hyperbolic case, 
i, e = — 1 — elliptic case. 



In particular, for stationary axisymmetric fields £ = z + 
ip, rj = z~ip, whereas for plane waves or for cosmological 
solutions £ = —x + t, n = —x — t, or these may have more 
complicate expressions in terms of x , x 2 . 



Dynamical equations 

The symmetry reduced dynamical equations for the 
action ([T]) can be presented in the form of real matrix 
Ernst-like equations for the string frame matrix variables 
- a symmetric dx <i-matrix Q, antisymmetric dx d-matrix 
B, a rectangular d x n-matrix A, the scalars $ and a: 

S = e 2 *||G ao ||, B=||B o6 ||, -A=|K (p) || 

which should satisfy the system of equations 

' if"dp(ad v £) -arf v {d lL £ - 2d n AA T )G- 1 d v £ = 0, 

< rfd^aduA) -arj^(d n £ ~ 2d tl AA T )g- 1 d u A = 0, 

, rTdn&a = 0, 

(5) 



where T means a matrix transposition and 

£=G + B + AA T , det^ = ea 2 e 2$ . (6) 

The equations ([3]) imply the existence of antisymmetric 
dx d potential B and dx n potential A defined as 



d^B = -eae. v ?Q-\d v B - d v AA T + Ad u A T )G~ x , 
duA = -ea£n' / g- 1 d u A + BduA. 



(7) 



The remaining (non-dynamical) part of field equations 
determines the conformal factor / in quadratures , pro- 
vided the solution of dynamical equations is found [l6| ■ 



Equivalent spectral problem 

As it was described in [l6j], the dynamical equations 
([5]) admit an equivalent reformulation in terms of the 
spectral problem for the four (2d + n) x (2d + n)-matrix 
functions depending on two real (in a hyperbolic case) or 
two complex conjugated (in the elliptic case) coordinates 
£ and rj and a free complex (" spectral" ) parameter w *E C 



(8) 



which should satisfy the following linear system for \& 
with algebraic constraints on its matrix coefficients 



2(w -£)<%* = U(£,7y)* 
2(^-7?)^* = V(£,r?)* 



U • U = U, trU = d 



V V 



V, trV = <2 

(9) 

It is necessary also that the system ^ will possess a 
symmetric matrix integral W (u>) such that 



* T W* = W (w) 
W^(w)=W (w) 



dW 

dw 



= n, n = 




where Id is dx d unit matrix and Q is (2d + n) x (2d + n)- 
matrix. We should impose also the reality conditions 



*(£,77,WJ) = *?,«>), W Q (W) = W H, W (3)(3) = In 

where W7 3 v 3 ) is the lower right n x n-block of W which, 
in accordance with (f5j)— (TTUj) . is a constant matrix and 
therefore the last condition in ([TTj) is pure gauge. 



Field variables and potentials 

As it can be shown by direct calculations (the detail 
will be published elsewhere), the conditions (|51)- (fTU|) im- 
ply, in particular, that W possess a special structure 

W = (w- j3)n + G, where 

fe^g- 1 -BQB + AA T BG + AA T A\ 



G = 



\ 



—GB + AA T 
A T 



G + AA T 
A T 



A 

In) 



(12) 
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and a — (£ — r))/2j, /? = (£ + r))/2; d x <i-matrix blocks 
Q and B are symmetric and antisymmetric respectively 
and, together with d x n-matrices A and A, these satisfy 
© - 0. (The corresponding expressions for U and V 
also have correct forms described in [jj].) This allows 
to calculate all field components and potentials for any 
solution (jSJ of our spectral problem. 



where fi* oJ Q no are some local regions on the complex 
plane near £ and r\ Q respectively such that the matrices 
U and V are holomorphic functions of (£, rj) in il^ o x £l Vo - 
Besides that, without any loss of generality we impose 
on the field components and auxiliary matrix functions a 
set of normalization conditions which provide unambigu- 
ous correspondence between local solutions of (JSJ and the 
solutions of our spectral problem. For this we set 



Global symmetries 

The spectral problem (|5j)- (fTTj) admits global symmetry 
transformations (including the discrete ones) 



U -> AUA- 1 , 
*^A*, W 



V -> AVA- 1 
■> (A T )- 1 WA- 1 



a 1 jia = n 



w 



(3) (3) 



In 

where the real constant matrix A is determined by two 
invariance conditions shown just above on the right. 
Some of these symmetries are not pure gauge and gen- 
erate physically different solutions from a given one. We 
use further two subgroups of pure gauge transformations: 




and A 



a a 









where u is antisymmetric d x d-matrix, a and a are dxn- 
matrices, Ld £ GL(d,M.), S n £ 0(d). These generate the 
following changes of the field variables and potentials: 



B' = B + u 
A' = A + a 
A' = A + ujA- 



and 



& = L^QLd, 

A! = LjAS n (15j 
A' = L- x AS n 



Another pure gauge transformation changes only the 
choice of the fundamental solution T], w) of © using 
arbitrary nondegenerate matrix C(w) 



(16) 



We need also a set of non-gauge discrete transformations 
which interchange in (fT3|) the fcth with (k + d)th rows 
and fcth with (fc + d)th columns of a matrix (1 < k < d): 



J (fc,d+fc): 



J (fc,d+fe) 



A (fc,<2+*0 



7 -l 



(17) 



*(£o,?7 , w) 



G„ = 








i, w H = (w- /3 )n + g 

G(to,rio) = Go, , 10 . 
B(£ o ,r, o ) = 0, (18) 

Mo,vo) = o, 

■A(£o,*?o) = 




where a = (£ Q — r/ a )/2j, f3 a — (^ +rj )/2 and dx d- matrix 



Go 



diag{e, 1, 
diag{e, 1, . 



■,1} - 
.,±1} - 



for a Q 7^ 
for a n = 



All these normalization conditions can be achieved using 
the transformations (|13[) - (fT6)) . Note, however, that in the 
case a = (as, e.g. in stationary axisymmetric case with 
the initial point chosen on the axis) , the finite expression 
for G D arises from the case a ^ as the limit for a — > 0. 
To get the non-degenerate Q in this case, one needs to 
use the transformation (fT5)l , (fTTf which may lead to a 
change of signature of Q Ql and after a construction of the 
solution, we have to perform the inverse transformation. 



General analytic structure of on w-plane. 

Global structure of on the spectral plane. As it is 
easy to see from the structure of © and (|18p , the normal- 
ized fundamental solution r], w) of the linear system 
([9]) possess in general only four singular points w — £ Q , 
w = £, w — r/o, w — rj which are the branching points of 
the orders (h, — 5, |, — 5) respectively for >&(£, 77, iu) and 
opposite for its inverse. To select a holomorphic branch 
of \P(£, 77, w), we make on the iu-plane two local cuts, 
one of which L + goes from w = £ to w = £ and another 
one L_ from w — rj Q to it) = rj. We assume that each 
of these cuts belongs to the corresponding local region 
S7+ = {w \ w = £ £ S\ o } and f2_ = {u> | w = rj £ ^l Vo }- 



The space of normalized local solutions 

The main object of our present consideration is the 
space of all (normalized) local solutions of ([5]), i.e. all 
solutions which components and potentials are holomor- 
phic (as functions of two complex variables £ and rj) in 
some local domains near a chosen "initial" point (£ Q , r) ). 
Therefore, we have to consider a general solution of (JHJ- 
([TTj) locally in £ and rj, but on the entire spectral plane: 




10 £ 



Fig.l Cuts on the w-plane in hyperbolic (e = 1) case. In the 
elliptic (e = — 1) case, the cuts are located symmetri- 
cally to each other with respect to real axis of ui-plane. 

Local structure of \l> near the cuts. It can be shown 
(see [l9| ) that the normalized matrix function \I/ and its 
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inverse possess near the cuts L± the following structures 

9 = f A; 1 t|> + (f , r/, w) • k+(w) + M+(£, to), u> G 0+ 
\ Al 1 tp_(^,?7,u;) • k_(io) + M_(£,7/,uj), to G 

x f A+ l+(w) • <p +(£, rj, w) + N+(£, T), to), w G 0+ 
I A- 1_( 

where "•" means mutrix multiplication and the functions 
A± are defined by the expressions: X±(w = oo) = 1 and 

A+ = V( w ~0/( w - lo)) A - = V( w ~ v)/(w - %); 

all other "fragments" of the local structures (fl9l) — d x 
(2d + n)-matrix functions k±(w) and tp ± ?7, w) , (2d + 
n) x d-matrix functions l±(w) and ip 77, u>) as well 
as (2d + n) x (2d + n)-matrix functions M±(£, 77, w) and 
N±(£,?7, iu) are holomorphic in i7 + or f2_ respectively 
and should satisfy there the algebraic relations 

k-N = 0, N-ip=0, l(w) = E^W-^w) -k T (w), 

M1 = 0, (p-M = 0, E» =(10 -£,)(«; -»/ ). 

(20) 

Here and somewhere below, we use instead of a pair of 
functions with indices "+" and "— " defined and holomor- 
phic in f2+ and f2_ respectively, a one function defined 
in + U 0_ and represented in Sl + and J7_ by the corre- 
sponding functions with the indices "+" or "— " . 



Monodromy data 



means that not all components of k±(i/j) are impor- 
tant in (j2"Tj) , and to reduce this ambiguity we should 
consider \t±(w) as functions which take the values in 
Grassmann manifold &d,2d+n(C). The last restriction 
on k±(iy) arises from the reality condition (jlip and it 
takes the form k±(u7) = \c±(w) in the hyperbolic case 
(e = 1) or k±(uJ) = k T (w) in the elliptic case (e = — 1). 
Just these Gd,2d+n-valued functions k.±(w) satisfying the 
above reality condition, parametrize unambiguously the 
whole space of local solutions of ([5]) . 

To specify the monodromy data as arbitary inde- 
pendent functional parameters, we chose the affine 
parametrization of the elements of Gd,2d+n(C). For 
this, we choose c±(w), so that any d linear independent 
columns of k + (w) as well as of k_ (w) will be reduced to 
a unite matrix Id . If these d linear independent columns 
can be chosen to be located at the same positions in 
k + (w) and k_ (w) (what is the generic situation, though 
there is some exceptional subspace of solutions for which 
k±(u>) can not be presented in the form given below, 
and these solutions should be considered separately) , the 
transformations (IT7l) can be used to locate these columns 
at first d positions. Then k±(w) take the forms 

k±H = u±H, v±(«;)} (22) 

where two d x d-matrix functions u±(ui) and two d x n- 
matrix functions v±(w), which are holomorphic in J7 + 
and 0_ respectively and satisfying the reality conditions, 
play the role of "coordinates" in the infinite-dimensional 
space of local solutions for (symmetry reduced) bosonic 
sector of the effective heterotic string theory. 



The coordinate-independent functions k±(iy) and 
l±(w) play a very important role in our construction. 
These functions determine the monodromy of >& on the 
cuts L± because the analytical continuation of along 
a simple path which goes from one edge of a cut L + or 
L_ to its another edge, leads to the linear transforma- 
tions of ^ — > • T± (w) with the matrices T± (w) which 
structure follows immediately from (1191) : 



T± (w) =1—2 l±(io) • (k ± (w) • 

i±H) 1 -M«0 (2i) 

where one should choose only upper or only lower signs. 
Note also that T^_(w) = I. 

In what follows, we are going to characterize unam- 
biguously the solutions of by a complete set of in- 
dependent functional parameters in (|21[) which we call 
the monodromy data of a solution. For this we note 
that in view of (|2"0")) , l±(w) can be expressed in terms of 
k±(w) and therefore, T±(w) are completely determined 
by k±(w). On the other hand, it is easy to see, that 
remains unchanged after the transformations 

k±(u>) — > c±(w) ■ k±(iu), ip ±{ w ) A> ±( w ) ■ c ± 1 ( w ) 

where c±(w) are arbitrary non-degenerate d x d-matrix 
functions holomorphic in f2 + and f2_ respectively. This 



Linear singular integral equations 

In accordance with the well known theorems of com- 
plex analysis, the analytical structure of 77, w) allows 
to present it and its inverse as the Cauchy-type integrals 



1 



i 



* = I + — / dC, * 1 = I + — / l - ± d( 



c 



c 



(23) 

where [V&J^ and [>& denote the jumps (i.e. half of the 
difference of left and right limits) of and 1 at the 
points (ei. The "continuous parts" {*}c and {\E' _1 }c 
(i.e. half of the sums of left and right limits) of '3/ and 
'3/ 1 on L are determined by the integrals of the form (j23|) 
in which, however, we put w = r G L and the singular 
integrals should be considered as the Cauchy principal 
value integrals (Sokhotski-Plemelj formula) [20|. From 
(fT9")l we obtain at the point £ G L = L + U L_: 

[*]« = [A _1 k ip (£, V, • k(C), {*}< = M(€, 77, C) 
[v-^c = [\] c 1(C) • v (f, v, 0, = N(C, v, 

(24) 

where we took into account that {X±} = on L±. Using 
these expressions for M and N in terms of Couchy prin- 
cipal value integrals in the relations (|20p k • N = and 
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M • 1 = 0, we obtain two systems 

1 f [A] c (k(r)-1(Q) 
-— f i 7 • <P (£ , V, d( = k(r) (25) 

L 

-lL,;,,, C ). HiM < = 1(r) (26) 

L 

which are the linear singular integral equations for matri- 
ces cp and \p , provided the monodromy data are given. 

The matrix equations (|25p and (|26[) are equivalent to 
each other, and for constructing solutions it is enough 
to solve only one of them. For pure technical reasons 
we choose ([2"5j) as the basic one. The theory of such 
equations is well developed [2(|, and it allows to show 
[19| that the solution of this integral equation (at least, 
for (£, rj) close enough to (£ ,f7o)) always exists for any 
choice of the monodromy data {u(u>), v(iu)}. 

Calculation of the field components 



Cauchyprincipal value integrals (Poincare-Bertrand for- 
mula [20] ) which takes in our context the form 

-/-I^MCR = B(r) & A(r) = Lf^k B (QdC 

ITT J Q — T ITT J Q — T 

L L 

for any function A(t) which satisfy the Holder condition 
on L and is bounded on its ends, as well as the elementary 
theory of residues which allow to reduce these integral 
equations to an algebraic system. 

Many physically important stationary axisymmetric 
solutions (such as, e.g. black holes in 4D case or black 
holes, black rings, etc. in 5D case) belong to this class 
and the integral equation method presented above allows 
to construct their multparametric generalizations. 

We note also that the fields with analytically not 
matched monodromy data (u + (u>) ^ u_(u>) and/or 
V-\-(w) ^ \—(w)) also include physically important types 
of fields, e.g. colliding waves or cosmological models 
which possess the singularities at a = 0. Methods for 
constructing such solutions also can be developed Q. 



For any given monodromy data and the corresponding 
solution (p (£, r), w) of the linear integral equation (|25|) 
the field components can be determined in quadratures. 
In accordance with (1231) and (l24l). for w — > oo we have 



* 1 = I -w X R 



R 



1 



Mc i(0-<p(£,r?,CHC 



where (2d+n) x (2d + rt)-matrix R(£, rj) is determined by 
the integral over L + U L_ . Then we obtain from ([8l)-(fTT]) 

U = 2a ? R, V = 2<9 J( R, W = W (w)-fJ-R-R T -rj 

Due to (fT2|) , all field components can be calculated alge- 
braically in terms of the components of the matrix R. 



A class of analytically matched monodromy data 



Solutions with analytically matched monodromy data 
u + (w) = u_(to) = u(w), v + (w) = V-(w) = y(w) (27) 

possess the following important properties: 

(i) the condition (|27|i is equivalent to regularity of the 
boundary a — (e.g., of the "axis of symmetry"); 

(ii) if we locate the initial point (£ OJ r/ ) on some part 
of the "axis" a = 0, the analytically matched mon- 
odromy data can be expressed in terms of the values of 
the Ernst potentials on this part of the axis (provided 
Go = diag{e, 1, . . . , ±1} - see the comments to ([IB])): 



{P) 2(0 ' nP) 08 -AO 



(28) 



(iii) for any choice of u(w) and \(w) as rational func- 
tions of w the corresponding solution of the field equa- 
tions can be found explicitly using the inversion rule for 



5D Minimal supergravity 

The equations for the bosonic fields in 5D minimal 
supergravity which are determined by the action 



F A F A A 



coincide with D = 5 and n — 1 bosonic equations of 
heterotic string gravity (fTJ) if we impose there the con- 
straints $ = and Habc = * Fabc with a subsequent 
rescaling Fab — > (2/V3)Fab- In our present context, 
these constraints are equivalent to the relations 

det£ = ea 2 , B ab = -e abc A c , B ab = -e abc A c (29) 



where e abc is the Levi-Civita symbol on the orbits of isom- 
etry group with coordinates {x 3 , x 4 , x 5 }; e 345 = £345 = 1. 
Thus, the space of solutions of 5D minimal supergrav- 
ity is embedded into the space of solutions of heterotic 
string gravity. Therefore, to construct the solutions for 
5D minimal supergravity, we have to impose the appro- 
priate constraints on the choice of monodromy data. This 
can be done directly for the class of analytically matched 
data taking into account (|2"5|) . however, the general class 
of monodromy data needs more complicate analysis. 



Vacuum gravity in D dimensions 

More strong restrictions imposed on the general mon- 
odromy data for heterotic string gravity in D dimensions 

detg = ea 2 , 6 = 0, ^ = (30) 

leads to the space of vacuum solutions in D dimensions. 
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Solutions and their monodromy data 



Concluding remarks 



We consider here some very simple examples of 5D 
solutions. The first of them is the Minkowski metric 



ds 2 



-df + dpi + dp 2 2 + pfdip 2 + p\dvb 2 



(31) 



with 



PXP2, P = Z 



Pi 



normalization on the "axis" pi 
(analytically matched) monodromy data for ([3T 



With the point of 
at pi = \J —2z Q , the 



are: 



diag{0, -—, —}, 







The four-parameteric solution for a charged rotating 
black hole in 5D minimal supergravity [21] also possess 
analytically matched monodromy data with one pole: 
u(w) = u + Ui/ (w — hi) and v(w) = vi/(w — hi) where 
the matrices Uo and Ui , the vector Vi and the scalar h are 
real and constant. For simplicity, we show their structure 
for non-rotating case - 5D Reissner-Nordstrom solution: 



u(w) = ( 





v(w) = ( 




where the parameters do, bo, Co, do, po and h depend on 
the mass m and charge q(— s m) parameters of a black 
hole and the position z — zq of the point of normalization: 



a 



m[m-2z (l + 2s 2 )] 



bo = -- 



4zg 
1 



d = 



2(30-/1) 
h = m(l + 2s 2 )/2 



Co 



Po 



(2z — m) 2 — 8ms 2 Zo 
2(z - h){Az 2 -m 2 ) 

,2 



Po = 



For a conclusion we note, that the infinite hierarchies 
of multiparametric families of solutions which can be con- 
structed explicitly using the approach described above for 
pure vacuum gravity or heterotic string gravity in D di- 
mensions as well as for 51? minimal supergravity include 
many known solutions and allow to construct their mul- 
tiparametric generalizations. In particular, this extends 
considerably the classes of soliton solutions in these grav- 
ity models. However, a lot of work remains and is neces- 
sary to identify the most physically interesting solutions 
and to analyse their physical and geometrical properties. 
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